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We describe a numerical model for the interaction of light with large raindrops using realistic 
nonspherical drop shapes. We apply geometrical optics and a Monte Carlo technique to perform 
ray traces through the drops. We solve the problem of diffraction independently by approximating 
the drops with areaequivalent ellipsoids. Scattering patterns are obtained for different polarizations 
of the incident light. They exhibit varying degrees of asymmetry and depolarization that can be 
linked to the distortion and thus the size of the drops. The model is extended to give a simplified 
long-path integration. © 2008 Optical Society of America 



1. Introduction 

Large raindrops arc deformed from spheres because of 
aerodynamic pressure as they fall at terminal veloc- 
ity. 1-3 Beard and Chuang's 4 (and Chuang and 
Beard's 5 ) model results for the flattening at the base 
of the drops are used in this study. Drop shapes are 
specified by 



r{6) = a[l 



10 



j cos(n0)] 



(1) 



where a is the radius of the undistorted sphere, lo- 
cated at the center of mass of the drop; 9 is the polar 
elevation with 8 = 0° pointing vertically downward 
(Fig. 1) and c„ arc shape co-efficients given in Ta- 
ble 1. Most previous research is restricted to spherical 
drops 6-8 or ellipsoid particles. 9 Mackc and Groklaus 10 
used the same distorted drop shapes as the present 
study for backscattering applicable to lidar measure- 
ments of rainfall intensities. In the current study we 
are interested in using drop distortion in a forward- 
scatter mode so as to measure drop size and, for long 
paths, rainfall intensity. A complete description of the 
model, including its application to a long path inte- 
gration, can be found in Ross. 11 
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2. The Ray Trace Method 

We used a Monte Carlo approach in which a large 
number of photons are traced through the drop using 
geometrical optics. Diffraction is accounted for ana- 
lytically. For our simulation we chose a wavelength of 
A = 650 nm, and thus the size parameter x = 2ira/\ 
is in the range 5-10 3 <a;<4-10 4 for the drop sizes 
in Table 1. Glantschnig and Chen 6 found that geo- 
metrical optics produced good results in comparison 
with the rigorous Mie theory for x > 20 and scatter- 
ing angles < 60°. Macke et al. 12 found this condition 
to be x > 60. In a different study, 13 computations 
were carried out over a wide range of size parameters. 
The results showed that the deviation from Mie the- 
ory is less than 1% for size parameters x s=s 10 and 
m = 1.33, confirming the validity of our geometrical 
optics approach. 

A. Forward Scattering 

Table 2, adapted from van de Hulst, 14 gives the forward- 
and backward-scattered intensities for spherical scat- 
terers and polarizations perpendicular and parallel to 
the scattering plane (polarization 1 and 2, respec- 
tively). Over 91% of light of polarization 1 and more 
than 97% of polarization 2 is forward scattered. As 
can be seen from Table 2, upward of 99.5% of the total 
forward-scattered light for both polarizations emerges 
from the first interface after simple reflection (p = 0) 
and from the second interface after twofold refraction 
(j> = I). The fraction of forward-scattered intensity in 
the p = and p = 1 rays increases for more distorted 
drops (see Section 4), and so it is sufficient to con- 
sider only the contributions of these rays to forward 
scattering. 
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Table 1: Shape co-efficients for cosine distortion fit (Eq. (1)) for drop radii between 0.5 and 4.5mm from Chuang and Beard. 5 
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B. The Procedure 

Fig. 2 shows the flow diagram for the ray trace model. 
The model is run with N = 10 8 photons. Runs for dif- 
ferent polarizations of the incident light use the same 
set of random photons to exclude statistical variations 
as a possible cause for differences in the scattering pat- 
terns. The photons propagate parallel to the x axis 
from negative x toward the drop, and a point P on 
the drop's surface can be represented in Cartesian co- 
ordinates by 



"(<?) 



— sin (p sin ( 
cos <fi sin 9 
— cos 9 



(2) 



where <f> is the azimuth angle, increasing clockwise 
from the y-axis. The surface normal at P is 



N 



dP dP 

~d9 X ~d6 



(3) 



Using Snell's law and the Fresnel equations, we can 
now trace the photons through the drop. The detec- 
tor coordinates (6,$) are defined similarly to (0,0), 
and arriving photons are resolved into angular bins of 



Table 2: Separation into forward and backward scattering for 
polarisations 1 (first number in each column) and 2. 
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1° resolution in and Each photon generated is 
assigned a weight of one. The intensity fractions for 
p = and p = 1 are stored in an intensity matrix I 
with O and $ as row and column indices, respectively. 

3. Diffraction 

Analytic solutions for diffraction from distorted drop 
shapes do not exist. Here, drop shapes are approxi- 
mated by area-equivalent ellipsoids, for which the diffrac- 
tion pattern is known. This approximation is veri- 
fied when the exact diffraction problem is numerically 
solved for one of the larger drops. Results showed 
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Fig. 2: Flowchart for the ray trace. 

negligible differences at the angular resolution used. 
Given the range of size parameters in this study, diffrac- 
tion scatters energy equal to refraction and reflection 
combined. 14, 15 i.e., the extinction efficiency is 2. 

A. Diffraction from Spheres 

The framework for the elliptical approximation is based 
on diffraction from a spherical aperture: 16, 17 



By integrating relation 5, we write the fraction of the 
total intensity contained within a cone of angle * as 

= 1 - ^ 2 (afcsin(*)) - ^(ak sin(*)) , (6) 

where ^fo is the Bessel function of order zero. 18 

If the energy flux through a ring element of area 
27rr 2 sin* d<J (Fig. 4) is L'(*) d*, then the fractional 
flux through a detector clement of angular dimensions 
(d6,d<I>) and area r 2 sin 6 d9 d$ is 



£a exp [i(cut — kr)] 



e =o 9=o 
x cos(6> — £)] g dg d9 , 



exp [i(kgq/r) 
(4) 



where <§ is the electric field amplitude, giving intensity 
/ oc |<?| 2 . £a denotes the source strength per unit area 
of the aperture; and r the distance between its center 
and the detector (Fig. 3). The solution is 



I (X 



' f x {kaq/r)~ 


2 r 


kaq/r 





(/casing) 



ka sin * 



(5) 



where a is the radius of the aperture, k is the wave 
number 2tt/X, and #i is the Bessel function of the 
first kind and order one. The solution does not de- 
pend on the angle £ due to the circular symmetry 




Fig. 3: Co-ordinate system used for the diffraction integral in 
Eq. (4). 



sin 6 d6 d$ 
27r sin * 



(7) 



which can be expressed in terms of detector coordi- 
nates because cos* = sin 6 sin $. After this substitu- 
tion and further simplifications, the final result is 



A9A# 



2 fl(aki) sin 2 9 sin $ 



d6 d$ (8) 



1 /2 

where £ stands for [cos 2 + cos 2 $ sin 2 ®Y<\ The in- 
tervals AO and A<I> are appropriately chosen to cover 
the size of the bin. By analogy with the intensity 
matrix introduced in Subsection 2.B, this yields the 
181 x 181 intensity matrix I for diffraction, which we 
find numerically. 

B. Approximating the Drop Shapes with Ellipsoids 

The Fraunhofer diffraction integral for the distorted 
drops has no analytic solution as g is a function of 9. 
Hence we approximate the drops with arcaequivalcnt 
ellipsoids. The cross-sectional area of the drops can 
be calculated from 



A= \i 







10 



2 



y] c n cos(n9) 



n=0 



1 + 2c + ci + 



10 r 2 

2 

n=l 



60 



(9) 



with the co-efficients c n from Table 1. The discrep- 
ancy between the approximation and the actual drop 
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To simplify Eq. (10), the following substitutions can 
be applied 
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Fig. 5: Comparison between the ellipses used (dash-dot line) 
and the actual drop shapes for drop diameters D = 2, 4 and 
6 mm. 



shapes is the largest for the bigger drops and vanishes 
as the size (and the flattening at the base) of the drop 
decreases. Fig. 5 shows some examples. 

Extending the circular aperture in one direction 
by a constant factor \i will cause the diffraction pat- 
tern to contract in that direction by the same factor. 
Because of the fi times larger area of the aperture, the 
intensity is /i 2 times the original intensity at each point 
mapped from the original pattern. 17, 19 Hence we can 
calculate the contribution to each bin from diffraction 
using elliptic obstacles from the results of a circular 
aperture. 

C. Verifying the Elliptic Approximation 

In this subsection the diffraction integral is evaluated 
numerically for one exact drop shape. The results are 
used to validate the elliptic approximation. The con- 
stant radius g in Eq. (4) is replaced by the cosine dis- 
tortion fit from Eq. (1), and the variable of integration 
is changed from g to a, the radius of the undistortcd 
sphere. Eq. (4) becomes 



<?((?, E) = 



£a exp [i(u>t — kr 



R lis 



x cxp 
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The integrand in Eq. (14) is resolved into its real and 
imaginary components 



*(<?,E)=C 



x£i?sin(x£#) + cos(x£R) - 1 



d6 



?=o 

J x 2 

9=0 



(15) 



which can be evaluated numerically for a given pair of 
detector coordinates (q,£). To return to the original 
detector coordinates (0, $) it can be shown that 

-cos(0 -S) = cos 9 cos9 - sin 9 cos $ . (16) 

r 

Making the replacements in Eq. (10) and (11) leads to 
the same integral as in Eq. (14) except that \ now rep- 
resents k times the right-hand side of Eq. (16) instead 
of the left-hand side. 

Because / cx |<?| 2 , only the real part of Eq. (15) 
needs to be solved. Fig. 6 shows the results obtained 
for a distorted drop with radius a — 3 mm. The 
locations of the first three minima exhibit a clear el- 
liptic symmetry. The dotted lines emerging radially 
from the center represent the lines along which the in- 
tegrations were carried out. The dotted circles with 
numbers next to them give the lines of constant \& as 
percentages of a degree. The resolution for the integra- 
tion was 10,000 points/deg, which gives approximately 
500 points from the center of Fig. 6 to the outermost 
circle. 

A continuous plot of the field amplitudes <£(q, 0°) 
and (o(q, 90°) is shown in Fig. 7. The eccentricity 
of the pattern is obvious. The first minimum in the 
S = 90° direction is considerably closer to the center, 
indicating the contraction in the horizontal. A simi- 
lar comparison between ${q, 0°) and ${q, 180°) (not 
shown) yields two identical curves within the margin 
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Fig. 6: Location of the first three minima in the diffraction 
pattern from a drop with radius a = 3 mm. 

of error, confirming the symmetry of the pattern about 
the horizontal. An interesting feature is the somewhat 
elongated sixth maximum in the £ = 0° direction. 
This feature returns in a similar shape with every 11th 
maximum. 

A more quantitative account of the possible flat- 
tening in the pattern for the 3-mm drop can be ob- 
tained from X(q) = (g(q, 270°) - 90°))/<?{q, 90°), 
which gives a measure for the relative deviation from 
horizontal symmetry. Values do not exceed 0.4% and 
are only notably over 0.1% within 0.1° of the central 
forward direction. We believe that this fully justifies 
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Fig. 7: Logarithmic intensity distribution for the directions (a) 
S = 90° and (b) E = 0°. 



the elliptic approximation, especially given the applied 
angular resolution of 1° and also considering that the 
distortion decreases for smaller drops. The same com- 
putations can be carried out along the £ = and 
£ = 180° directions. Values in the absolute deviation 
do not exceed 4 • 10~ 6 , which is the region of error 
applied during the numerical integration and can thus 
be neglected. 

4. Results 

The results from the computer model are presented 
separately for the ray trace (no diffraction) and for 
scattering including diffraction. The scattering pat- 
terns are displayed by a sinusoidal projection method 
(which means that the meridians appear as sine func- 
tions), hence the scale is preserved only on the central 
meridian and along the lines of latitude. This projec- 
tion method also renders the region around the central 
forward direction without distortions. 

A. Scattering from Single Drops - excluding Diffraction 

The phase is omitted because the path lengths of the 
photons arriving at a particular bin can differ by sev- 
eral tens of wavelengths. These effects cancel out for 
a large number of photons and do not contribute any 
valuable information to the pattern. 

Fig. 8 shows the scattering patterns for a 3-mm 
drop and 10 8 photons. In Fig. 8(a) the incoming light 
was unpolarized. The pattern exhibits a symmetry 
about the central meridian that is to be expected from 
the rotational symmetry of the drops about their ver- 
tical axis. This can also be observed in the corre- 
sponding polar plot of Fig. 9(a) that shows the inten- 
sity distribution along the = 82° line that contains 
the absolute maximum. The concentric semicircles are 
lines of constant intensity, and the numbers give the 
corresponding intensities as powers of ten. Fig. 9(b) 
shows the distribution along the $ = 90° meridian. A 
distinguishing feature is the lack of symmetry about 
the equator. The absolute maximum has shifted from 
the equator (0 = 90°) to = 82°. It is accompanied 
by another local maximum at = 130°. For smaller 
drop sizes, the absolute maximum shifts closer to the 
(90°,90°) direction again. 

Another prominent feature in Fig. 8(a) is the 
sharp gradient, creating a steep intensocline that sep- 
arates the central area of relatively high intensity from 
the outer, low-intensity regions. It is identified by 
the narrowly spaced contour lines. The model was 
used to record the intensities from the rays p = and 
p = 1 separately The light that arrives outside the 
intensocline has mostly undergone a single reflection 
at the first interface, whereas the high-intensity area 
is mainly lit by the twice-refracted rays. 
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Fig. 8: Contour plots of the intensity distribution obtained from a 3 mm drop, for: (a) unpolarised light, (b) y- and (c) z-polariscd 
light. A plot of the difference in the obtained intensities for y-polarised and 2-polarised light (non-logarithmic) is shown in (d). 
Note that the gray scale below (c) applies for (a) to (c) while (d) being the only non-logarithmic plot has a scale of its own. 



Fig. 8(b) shows the pattern for the same drop but 
with incident light polarized parallel to the y axis. 
There are two distinct minima near the equator at 
extreme angles for $. This can also be seen on the 
corresponding polar plot Fig. 9(c). The sum of contri- 
butions from all angular bins is slightly (« 1%) larger 
than for unpolarized light. The opposite is true if 
the incident light is polarized parallel to the z axis as 
shown in Fig. 8(c): The distinct minima arc on the 
$ = 90° meridian, symmetrically at 9 = 90° ± 74°. 
Light that is reflected into the angular bins containing 
the observed minima in Figs. 9(c) and 9(d) originates 
from points on the drop's surface where the angle of 



incidence is at the Brewster angle of a p = 53.06° and 
where the plane of polarization of the incident light is 
parallel to the scattering plane. 

Fig. 8(d) shows the difference between the re- 
sults obtained for y-polarized and z-polarized incident 
light. This plot shows the absolute deviation, and the 
minima that can be seen in Figs. 8(b) and 8(c) do 
not appear. There is, however, a significant differ- 
ence in the high-intensity area, inside the intcnsocline 
and especially near the absolute maximum. Note that 
the scale of this last plot is nonlogarithmic to allow 
for both positive and negative differences. A high- 
intensity ridge that extends from 6 = 78° to 9 = 86° 
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Fig. 9: Polar plots for a 3 mm drop showing the intensity distribution along: (a) 8 = 82° and (b) <E> = 90° , both for unpolarised 
light; (c) © = 82° for j/-polarised light and (d) <E> = 90° for z-polarised light. The graphs in (e) and (f) show the difference 
between the results for y- and z-polarised light along = 81° and <E> = 90° respectively. Negative values are plotted into the 
lower half of the circles. 



on the central meridian contains values that exceed 
the maximum value covered by the scale (compare 
with Fig. 9(f)). The actual maximum in Fig. 8(d) 
is 9200 and is located at 9 = 81°. The maximum 
difference of 9200 photons at the central maximum is 
equivalent to 12.5% in relative terms, which explains 
why this feature does not appear on the individual 
contour plots in Figs. 8(b) and 8(c). If the pattern in 
Fig. 8(d) is divided by the (nonlogarithmic) intensity 
from Fig. 8(a) (result not shown), the extrema outside 
the intensocline visible in Figs. 8(b) and 8(c) can be 



observed. The intensity difference at these extrema is 
nearly twice the intensity value that is obtained there 
for unpolarized incident light. If the incident light is 
unpolarized, the reflectance R at a p is half of the value 
of R±. Because R» is equal to zero, the difference be- 
tween the y- and z-polarized intensity is equal to R± . 
Dividing R± by R for unpolarized light thus always 
yields 2. The finite size of the angular bins allows in- 
clusion of rays for which the angle of incidence differs 
slightly from a p . This should always keep the abso- 
lute value of the extremes slightly below 2. A drop 
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Fig. 10: (a) The fraction of the total incident intensity that is forward scattered and (b) the fraction of all incident rays that arc 
entirely internally reflected (regardless of their intensity). 



of radius a = 0.5 mm, for example, has its extreme 
values at 1.98 and -1.95, and the corresponding val- 
ues for spherical drops are ±1.95. The reason why 
this is not observed for larger drops along the cen- 
tral meridian might be a combination of the varying 
curvature along the vertical axis and statistical fluctu- 
ations because the bins at such extreme angles contain 
only approximately 35 photons out of the total of 100 
million. 

Fig. 9(e) gives a clearer image of what happens 
along the 9 = 81° latitude. The scale is again log- 
arithmic to accommodate the large spread of values. 
The pattern is dominated by the central maximum 
and two broad maxima on each side. For smaller drop 
sizes the central maximum becomes less dominant and 
is almost nonexistent for a 1-mm drop. The two pe- 
ripheral maxima, however, remain almost constant in 
size and location. 

Fig. 9(f) shows the intensity distribution along 
the central meridian. It is characterized by closely 
spaced maxima and minima. There is a general de- 
crease in the number of significant extrema to be ob- 
served with decreasing drop size, accompanied by a 
general increase in overall symmetry 

Over 99.5% of the forward-scattered light from 
spherical scatterers emerges from either single reflec- 
tion (ray p = 0) or twofold refraction (ja = 1). The 
assumption was made that this would be similar for 
distorted drops, thus justifying rays with p > 1 being 
neglected. 

Fig. 10(a) shows the percentages of forward-scat- 
tered intensity in p = and p = 1 for unpolarized inci- 
dent light and drop sizes ranging from a = 0.5 mm to 
a = 4.5 mm. The fraction of forward-scattered inten- 
sity increases with drop size, reaching a maximum for 
drops with a = 3.5 mm. The rapid decline for larger 



drops can be explained from Fig. 10(b), which shows 
the percentage of photons that are internally reflected 
and thus discarded by the model. It can be seen that 
internal reflection can be neglected for most drop sizes 
and becomes noticeable only for large drops. 

To be able to give a quantitative comparison with 
some tabulated values, the model was used to find the 
scattering pattern that would be produced by unpolar- 
ized light that is incident on a sphere of radius 0.1 mm. 
The result is given in Table 3. It compares the fraction 
of energy that is scattered forward in rays p = and 
p = 1 from the computer model and Table 2. These 
differ by less than 0.2%, justifying the model and its 
assumptions. 

B. Scattering from Single Drops - including Diffraction 

The diffraction pattern is dominated by a narrow diffrac- 
tion peak in the central forward direction. The width 
of this diffraction peak depends on the size of the scat- 
tcrcr. The effects of diffraction are therefore demon- 
strated for two extreme drop sizes: a 0.1-mm spherical 
drop and a 3- mm drop. Figures 11(a) and 11(b) show 
the combination of diffraction and ray trace together 
with the results from the ray trace alone. The cen- 
tral maximum for the 3-mm drop is much sharper and 
slightly higher than for the smaller drop. The solid 
and dashed curves only differ noticeably near the cen- 
tral diffraction peak, whereas the intensity at other 
angles is almost unchanged. This is significant for 



Table 3: Comparing results from the Computer Model and van 
de Hulst. 14 



Source: 


model 


van de Hulst 


% forward scattered 
in p=0 and p=l 


93.73 


93.89 
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Fig. 11: Superposition of the polar plots of the logarithmic intensity including all components (continuous line) and without 
diffraction (dashed line) along the 90° latitude for (a) a = 0.1mm and (b) a = 3 mm. 



the design of an instrument to measure drop sizes. If 
the detectors are placed well off the central diffraction 
peak, the results can be simulated with contributions 
from the ray trace only. This is illustrated in a more 
quantitative form in Fig. 12. It shows the percentile 
change of the intensity distribution from the ray trace 
after diffraction is added. For the large 3-mm drop, 
the pattern is affected only by diffraction for angles 
of \& < 1° if we accept a 1% margin of error. For 
the smaller drop the same criterion would yield the 
condition * < 10°. 

5. Practical Application 

If a large number of individual field measurements is 
averaged, the underlying pattern linking the rainfall 
rate with drop diameter shows an exponential charac- 
ter as first proposed by Marshall and Palmer. 20 They 
suggested the following empirical relationship, the so- 
called MarshallPalmcr distribution: 

N(D) = N a exp (—AD), N = 0.08 cm" 4 ; (17) 

where N(D)dD is the number of drops per unit vol- 
ume having diameters between D and D + AD. The 
intercept parameter Nq gives the intercept with the 
ordinate for D = 0. The distribution depends entirely 
on the parameter A that is determined by the rainfall 
rate R: A = 41i?~ 021 where R is in millimeters per 
hour and A is in inverse centimeters. 

We can now establish a function P(D) to repre- 
sent the cumulative probability of finding a drop with 
a diameter between and D in a large enough sample: 



ization), one obtains the following area distribution: 



P(D) = 1 - exp (—AD) 



(18) 



If the drops are approximated as spheres of area ir/AD 2 
(which is possible because the number of drops de- 
creases rapidly as the distortion increases - see Eq. (18)) 
and by integrating Eq. (17) (while neglecting the con- 
stants Nq and ir/4 that both cancel out in the normal- 



ly 



J D' 2 cxp(-AL>')dD' = -exp(-AD) 



D 2 2D : 
X + A 2 " + A 3 / ' A 3 



(19) 



Normalisation yields: 

D 



J D exp (-AD')dD' 



o 



/ D 2 cxp(-AD)dD 



exp(-AD) 



,^! + D A + l 



(20) 



Fig. 13 shows the result for different rainfall intensi- 
ties. For R = lOmm/h, more than 50% of the in- 
cident light falls on drops with D > 1 mm, which is 
where drops start showing first signs of distortion. For 
heavier rain, more larger drops are present; and, as in 
the example of a R = lOOmm/h event, approximately 
80% of the incident light falls on distorted drops. This 
suggests that drops with D < 0.2 mm do not con- 
tribute significantly to the scattered intensity and can 
be neglected, which would produce only a 5% error 
even for light rain with predominantly small drops. 

A. Multiple Scattering 

The optical depth of a path of length z can be defined 



as: 



N e Cz, 



(21) 



where N e is the number of drops per unit volume in 
the sample and C is the extinction cross section. As 
was discussed in Section 3, the extinction cross section 
approaches a limiting value of 2 A (where A is the area 
of the scatterer) for particles with x 3> 1, and thus 



2ira z N e 



(22) 
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Fig. 12: The influence of diffraction on the overall scattering pattern. The percentile change between the ray trace only and the 
ray trace plus diffraction is shown along the 90° latitude for (a) a = 0.1mm and (b) a = 3 mm. 



where the drops have again been approximated as 
spheres. N e is given by 



oo 

N e = J N exp (—AD) dD . 



(23) 



Making the substitution, we find the area-distribution 
from Eq. (19) in a slightly different form: 



oo 

t z = 2irN z J a 2 exp (—AD) dD 
o 

oo 

= 4ttN z J a 2 exp (— 2Aa) da 



TTiVr, 

A 3 



(24) 



(25) 



(26) 



Fig. 14 shows the optical depth as a function of rain- 
fall intensity for four different path lengths. Van dc 



Hulst 14 argues that single scattering prevails if t z < 
0.1; corrections for double scattering are necessary for 
0.1 < t z < 0.3, and multiple scattering is observed 
for even larger values of t z . Other sources 21 maintain 
that multiple scattering can be neglected for t z values 
of up to one. 

B. Numerical Simulation 

The rainfall event is simulated when a large volume 
V is created above an imaginary laser beam and filled 
with drops according to the MarshallPalmer distribu- 
tion for a given rainfall intensity. Only drops in the 
size range 0.2 mm < D < 8 mm were considered for 
the reasons outlined after Eq. (20). The total number 
of drops in V is given by 



0.8 



N v = V No exp -AD dD 



VNp 
A 



0.02 



exp (-0.02A) 
(27) 
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Fig. 13: Area presented to the incident light with drop diameter 
for three different rainfall intensities. 
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Fig. 14: Optical depth as a function of rainfall intensity for four 
different path lengths. 
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Fig. 15: Dependence of average received light intensity on rain- 
fall rate R for (0,*) = (90°, 75°) and (©,<!>) = (90°, 27°). 



with the approximation exp (— 0.8A) sa yielding neg- 
ligible errors even for high rainfall intensities. 

It has been shown experimentally (e.g., Ref. 22 ) 
that, for a constant mean rainfall rate R, the number 
of drops is Poisson distributed around the mean. Ny 
is therefore used as a mean value to produce a Pois- 
son deviate Np. Once the volume is filled with Np 
drops, each drop is randomly assigned a size from an 
exponentially weighted distribution: 



D 



In a; with x G ]0, 1] 



(28) 



where x is a randomly generated number. The drops 
are then released and fall at their terminal velocities 
(with the values from Gunn and Kinzer 23 ) through 
V . Every t seconds a snapshot is taken, recording the 
number of drops present in the beam and their sizes. 

C. Results 

The scattering patterns are obtained when the pat- 
terns from individual drops arc combined according 
to their abundance in the beam. The patterns from 
the individual drops are normalized with the cross- 
sectional area to obtain a constant intensity per unit 
area. The result is a time series that gives the varia- 
tions of the received intensity for each angular bin. 

As expected, the results show that the higher 
the rainfall rate, the higher the mean scattered in- 
tensity over the sampling period. This is shown in 
Fig. 15 for two different angular bins. The mean de- 
tected intensity is thus a conclusive indication of the 
rainfall intensity R. The degree of asymmetry about 
the equator can also be used to infer the rainfall rate. 
Fig. 16 shows the difference between angular bins at 
(6,$) = (90 o -A6,90 o ) and (9,$) = (90°+A6, 90°). 
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Fig. 16: Relative deviations from vertical symmetry in the scat- 
tering data for different rainfall intensities and different angular 
bins. The deviation is given as a percentage of the total received 
intensity at (0, $) = (9O°+A0, 90°). 



As can be seen, the amount of light received in the bin 
below the equator is considerably higher than for the 
corresponding element symmetrically above the equa- 
tor. The difference increases as expected with rainfall 
intensity and varies depending on the distance to the 
equator. 

Although the light intensity is a stable indicator 
of the rainfall rate, the effects of wind on the orien- 
tation of the vertical axis of the drop can be severe 
on the above symmetry considerations. Other factors 
that are likely to occur in real situations and that are 
neglected in this model arc coalescence, drop breakup, 
and oscillations, as well as multiple-scattering events. 
Possible realizations of field experiments are given in 
Ref. 11 

6. Summary 

The distorted shape of raindrops at terminal velocity 
is often ignored and approximated with cither spheres 
or ellipsoids. In contrast, the numerical model devel- 
oped in this study determines the scattering pattern 
for the true drop shapes. This was achieved through 
a combination of geometrical optics and a statistical 
Monte Carlo technique. The ray trace yielded scatter- 
ing patterns for different drop sizes and polarizations 
of the incident light. The shape of the drops could 
be inferred from the varying degrees of depolarization 
of the scattered light and the asymmetries observed 
in the overall scattering behavior. No experimental 
data are available for light scattering from distorted 
drops. However, the results obtained for small (spher- 
ical) drops match tabulated values confirming the va- 
lidity of the model. 

We have treated the problem of diffraction by 
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successfully approximating the drops as ellipsoids. The 
error incurred has been found to be negligible even for 
the largest drops examined. Adjustments for the re- 
sults from the ray trace were necessary only near the 
central diffraction peak whereas diffraction can be ne- 
glected for the remaining pattern. 

The asymmetries in the scattering behavior can 
be used for rainfall measurements with a laser and 
several detectors as the drop size distribution depends 
on the rainfall intensity and thus influences the over- 
all symmetry. Further details including the experi- 
mental realization and application to high-resolution 
measurements of rainfall structures can be found in 
Ross. 11 

This model was developed as part of a Diplo- 
marbeit (equivalent to a Master of Science project) at 
the Fachbereich Physik of the Freie Universitat Berlin 
conducted in close collaboration with the Physics De- 
partment of the University of Auckland that gener- 
ously provided the resources for this research. 
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